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Abstract 

Correlation functions of QCD currents with quantum numbers of nucleon and A- 
isobar are considered at finite temperatures. Corrections of order to the correlators are 
calculated and interpreted in terms of thermal mass shifts using a QCD sum rules type of 
argument. In both cases the masses decrease with T. 
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The study of QCD current correlators at finite temperatures has been of increased 
interest recently (for an extensive review, see Ref.[||). By studying their T-dependence one 
may get an insight on the changes in the hadron spectrum as the temperature approaches 
the point Tc ~ 150 MeV at which the hadron phase is expected to undergo chiral and 
deconfining phase transitions. One may try to extend the method of operator product 
expansion (OPE) from T = to the case of T 7^ and use it to evaluate the finite 
temperature correlators. The T-dependence of non-perturbative condensates is due to 
interactions with the heat bath. Since at low T the heat bath is dominated by pions, it is 
possible to use soft pion methods to obtain it[§, ^, ^. However to leading order O(T^) the 
temperature corrections to correlators may be determined on the basis of current algebra 
and PC AC without actually using OPE|^, ^, |^. In this order a mixing of correlators in 
parity and/or isospin takes place||8|, ^. The important point is that this mixing arises 
from direct scattering of the heat bath pions on the currents and no thermal mass shift is 
generated in this order. (A general argument regarding the absence of mass shifts in this 
order was given in Ref.|Q). It has been shown |l^ that the T-dependence of the relevant 
condensates exactly matches the mixing of the correlators leaving no room for mass shifts. 

The condensates mentioned above are due to Lorentz scalar operators in the OPE. An 
important new feature of OPE at T ^ is that operators with non-zero Lorentz spin s 
contribute to correlation functions because their average over the heat bath is non-zero. 
They contribute to higher orders in T^ and are new non-perturbative parameters to be 
determined dl, ^, ^. They cannot be obtained using soft pion approach. 

In a recent paper |11| we have shown that order T^ corrections to the vector and axial 



correlators may be obtained using dispersion relations for the amplitudes of deep inelastic 
scattering on pions and discussed their relation to the lowest s = 2 condensate, the energy- 
momentum tensor. The results were interpreted in terms of p and ai thermal mass shifts 
which turned out to be negative. In the present paper we consider T^ corrections for 
correlators with quantum numbers of a nucleon and isobar. There is no similarity with 
deep inelastic scattering in this case, but we can use OPE with account of the s = 2 
operators to obtain the corrections responsible for the mass shifts. 
Let us consider the thermal correlation function defined as 

CN{q,T) = i J d^xe"i^{TMx),m})T (1) 

where 

(r{.M.,7(0)}). = g--<"l^'''f'-?i°»'"^^l"> , (2) 

Z^nX'^l^ 1^/ 
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Here H is the QCD Hamiltonian, the sum is over all states of the spectrum, and the 
nucleon (proton) current is taken to be|12| 



ry = KC77X)757/.d'=e'^'' (3) 

where C = 7270 is the charge conjugation matrix. It is assumed that is space-like, 
= —q^ > 0, and is much larger than a characteristic hadronic scale, ^ ^c^j 
where Rc is the confinement radius, R^^ ^ 0.5 GeV. 

To determine the correlator Cat at low T one has to calculate the matrix elements in 
Eq.(^) over pions. In the chiral limit the pion momenta are of order T, and the leading 
O(T^) corrections are obtained neglecting the pion momentum in the one-pion matrix 
elements. In this order using PCAC, current algebra and doing the integral over the 
thermal pion phase space 

/ ^^^i ? (-''(pm^(-)^^m-''(p)) (4) 

gives the exact result 0, [l^ 

CN{q,T) = (l-e)Civ((?,0) -^75^^(9,0)75 (5) 

where 



^ = rVl6F^ (6) 

and = 93 MeV is the pion decay constant. The T'^/F^ corrections are due to the 
contact interactions of pions with the currents (Fig.l, a and b). There is no need to use 
OPE to calculate Cat in this order. Indeed, the T-dependences of s = condensates which 
contribute to this order exactly satisfy Eg. (p|) pO|] . Thus, in order T'^ the finite temperature 
correlator C]\f{q,T) is expressed in terms of the T = correlators. This is similar to the 
mixing of the vector and axial correlators!^, but in the nucleon case both Cn and ^^Cnj^ 
are contributed by the nucleon and its parity partner 7V(1535). The choice of the nucleon 
current in Eq.(^) allows to avoid the compensation of their contributions to the chirality 
violating part of Cat 1 13]. It is evident from Eq.(|5|) that there can be no nucleon mass shift 
in order T"^. The absence of the nucleon mass shift in this order was also demonstrated in 
Ref . by an explicit calculation of the vr self-energy correction to the nucleon propagator 
at T ^ 0. The mass shift is expected to appear in the next order, T^. 

The corrections to Cat are due either to two-pion matrix elements with zero pion 
momenta, or to one-pion matrix elements with non-zero pion momentum [11]. The former 
are of order T'^/F^. They come from the contact interaction of pions with the currents 
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(a) 



(b) 




Fig. 1. Contact terms to two loops. Double lines are T = correlators, crossed dashed 
lines denote thermal pions. 



in two loops (Fig.l, c through h), can be calculated without OPE, and contribute to the 
change in the coupling of the nucleon to the current (see Eq.(|l3|) and ([l5| ) below). The 
latter are of order T^/Q^, they arise from the pion scattering on the intermediate state 
in the correlator and may be interpreted in terms of a thermal mass shift. They will also 
break the Lorentz covariance preserved in order as is seen from Eq.(|5|) in which Cn 
has the form 

CNiq\T)=Ciiq\T)q + C2iq^T) (7) 

where q = q^^^i- The general form of the correlator in the rest frame of the heat bath is 
given by 

CN{q,T) = Ci{q,T)^oq^ + C'2{q,T)-fq + C2{q,T) (8) 

where C( 7^ — Ci and all C's depend separately on and |q|. We will consider the case 
of q = 0, (7o 7^ 0, when we are left with the same number of independent structures as at 
T = 0. 

In terms of OPE, there are two s = 2 operators in the leading twist which provide 
the T'^ /Q^ corrections to the correlator. They are the quark and gluon energy-momentum 
tensors 



'^Iv = l^iqifiDuq + qiuDf,q), {q = u , d) 



normalized according to 



{<p)Ku + + 0%Hp)) = ^Pf^Pu (10) 

where (7r(p)|7r(p')) = {2tt)^ 2E 6^^\p — p'). The contributions of the quark and gluon 
tensors to Eq.(|lO|) are about the samep^, Q. However, the gluon tensor enters OPE for 
the correlator being multiplied by as^Q"^) and we will neglect its contribution. In what 
follows we will also omit the anomalous dimension of 0-^^. 

The calculation of the quark tensor contribution to the correlator can most easily be 
done in the coordinate representation by cutting quark lines and using the second term in 
the expansion of the quark operator, q{x) = q{0) + xxD\q{Q) + .... 

Cutting one quark line gives the chirality conserving part. In case of d-quark this 
results in matrix elements of the form {■K{p)\dD\d\T:{p)) which are proportional to the 
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(a) 



(b) 



Fig. 2. Leading terms in OPE related to the nucleon thermal mass shift. The quark leg 
with a cross corresponds to x\Dxq{0). 



pion momentum p and vanish under the integration in Eq.(^). Cutting a u-quark line 
(Fig. 2a) gives 

i^^T^{p{px)x'^ + x{px)'^) (11) 

where M2 is defined according to {tt{p)\0^^ + 0'^^^\7r{p)) = '&M2P^Pu and is equal to the 
second moment of the pion structure function M2 = F2{x)dx^A\. 

Integrating over the thermal pion phase space in the rest frame of the heat bath with 
the help of 

J (2vr)3 2|p| exp(lpl) - 1 " ^"^^^ + ^9,o9uo)^ (12) 

and going over to the momentum representation, we get the expression to order for 
the chirality conserving structure proportional to 70^0 

Ci(Q2,r) = Ci(Q^0)-^TMnQ2 (13) 

15 

(order ^ and contact terms cancel in this structure). 

Similarly, the chirality violating part of the T^/Q^ correction is obtained by cutting two 
quark lines. Here, cutting the two u-quark lines gives terms linear in the pion momentum 
which vanish under the integration over pions, while cutting the d-quark and a n-quark 
line (Fig. 26) gives terms quadratic in the pion momentum 

24M2(e>x)2 , ^ 

After integration over pions, this gives in the momentum representation 
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C2{Q',T)=(l-2C--e]C2iQ^0) (15) 



Here (qq) = -{24:0MeV)^ is the usual T = quark condensate [|. In Eq.(|l3|) and (|l^) we 
have taken into account order /F^ corrections in a way similar to the case of vector and 



axial correlators! 11]. The full answer was again obtained by taking into account the initial 
(finite) state interaction between two pions (Fig. 2, g and h). The leading term in OPE 
for C2 (Q^ , 0) is proportional to {qq) . A straightforward check shows that the .^^ term in 
Eq.(|l^) agrees with the knownQ correction to the quark condensate. It is easy to see 
that the last terms in Eqs.(p^) and (15) are indeed of order T^/Q^ compared to Ci((5^, 0) 



and C2((5^,0), respectively. 

Now, we would like to interpret the corrections to the correlators in Eqs.(|l3[) and 
(|T^ ) in terms of the nucleon thermal mass shift. The functions Ci and C2 may be presented 
in the form of dispersion integrals 



1 r'^lmCi{s,T)ds 
^ Jo s + Q^ 

where a sufficient number of subtractions is implied and a standard model (lowest reso- 
nance + continuum) for the spectral densities is usually assumed 



ImCiis, T) = 7TX^iT)6{s - m%{T)) + 9(3 - soiT))pis, T) (17) 
At T = the nucleon pole contribution to the correlator is 

Cr^(Q^O) = A^^^ (18) 
where Aa? = Ai(0) = A2(0) is the coupling of the nucleon to the current, 

(0|r/|iV) = Xnv (19) 

and the nucleon spinor v is normalized as vv = 2m n. At T 7^ the general form of the 
nucleon propagator is 

^ (20) 



7ogo(l + OAf) - mNil + bN) 



^In deriving Eq.(|lJ) we used the following approximation for the matrix element over pion 
{Ti{p)\u'y^DuUuu\n{p)) = {Ti{p)\u'yf_iDi,u\-K(p)){uu) . An estimate of the accuracy of this factorization may 
be obtained by inserting a pion intermediate state. Using the scale of OPE /i ~ 0.5 Gev as the cut-off in 
the integral over the momentum, we get that this contribution is suppressed as /Stv'^ ~ 0.3. 
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where a^r and bjy are T dependent corrections. Multiplying the numerator and the de- 
nominator in the above formula by 7o(?o(l + cln) + 'TiAr(l + bf^) we identify the thermal 
shift of the nucleon pole as 



6mN{T) = mN{bN{T) - aN{T)), (21) 

While SmN{T) here is of order T^, the functions aiv(T) and bi\i{T) contain terms 
ajyiT) = + 0{T'^), bN{T) = g\£_ + O(T^), where qa = 1-27 is the usual axial coupling 
of the nucleon. They are due to the contribution of the ttN self-energy graph at T 7^ 0. 
On the other hand, Eq.(|l3|) and (^), which are exact, do not contain any terms involving 
g\. This seeming paradox was resolved in Ref.[10| in which it was shown that the g\ 
terms indeed cancel provided the vrA^ — > N scattering is accounted for in the spectral 
density of the correlator at T 7^ 0. Taking this into consideration we obtain the nucleon 
contributions to Ci and C2 



Ci{Q\T) = Af(0)- 



Q-^ + m% + 5m%{T) 

c.{Q\T) = m (1 - 2c - '-e) "^"^I7"f^^/f7;/^^ (22) 

V 3 / Q-' + m% + dmjf{T) 

where aj\[{T) = 0{T^) is the nucleon wave-function renormalization due to irX self-energy 
corrections to the nucleon propagator at T 7^ 0, where X ^ N. To relate 6m]\}(T) and 
aAr(T) to the corrections obtained before, we match Eq.(p2|) against Eq.(|l3|) and ( |T5| ) 
keeping terms linear in 5m n and a^. To suppress the contribution of higher states, the 
Borel transformation 

1 '^2^n^ d 



BCm= lim rryW^r (23) 

/n=M^=const 

is usually performed in QCD sum rules. Though we neglect the changes in the continuum 
induced by temperature, the Borel transformation is still necessary to kill the subtraction 



terms proportional to in Eq.(13). Applying it, we get 



rriN ~ 15A^ ^ ^ 2Mh^ ) 

= l5A^ 

We see that the mass shift is negative. If varies between 0.7 and 2 MeV^ (the stability 
range of the T = sum rules for the nucleon), 5m ^ changes by 30%, while a is less 
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stable. Taking M = tun, = 0.12 Q and = 1.2 • lO^^ GeV''lT|], we find that 
the mass shift is rather small: dniN ~ —2 MeV at T = 80 MeV. We note that up till 
this temperature our estimate for 5m]\f agrees with the result of Ref.|^], obtained using 
experimental information on the vrA^ forward scattering amplitude. 

Taking into account higher twist s = 2 condensates, when they become available, one 
could make the above predictions more accurate. (Note, however, that is known with 
an accuracy of 50%[|l2|). Also, to go to higher T in this approach it is necessary to known 
higher spin condensates which contribute to next terms in the low T expansion of Jm^y. 

We have also considered the case of A-isobar. In that case the current in question 
isP 

V^ = {u'^C^y)u'e^''^ (25) 

The corresponding correlator C^j^{q) has a number of tensor structures already at T = 
and is contributed not only by A, but also by the resonance with = ^ . However, 
the structures g^uQ and g^i, correspond only to = | states. At T 7^ independent 
structures proliferate and we again simplify the situation by putting q = 0. The calculation 
is quite similar to the case of nucleon. The pion matrix elements which give rise to the 
T^/Q^ corrections proportional to g^y are 

- ^^-^{yP^lXu + PuX^)x{px) - g^j,ux{px) ) (26) 
for the chirality conserving and 

^^^^^^{(Pi^Xu +PuXf,){px) - gf,u{pxf) (27) 

for the chirality violating structure. Integrating over the thermal pion phase space and 
going to the momentum representation we get 

45 

and 



T^lng^ -g^uqolo (28) 



45 Q2 i^y; 

Matching these corrections against the contribution of A to the spectral density, we get 
after the Borel transformation 



*See [0 for a discussion of existing estimates of AI2 
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1 + ... T 
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The mass shift is again negative. Substituting = 2.5 • 10 ^ GeV^[12| and putting 
M = rriA, we see that mass shift is much smaller than the nucleon one: 5mA ~ —0.25 
MeV at r = 80 MeV. 

To summarize, we have used the only presently available non-scalar condensate of spin 
2, the energy-momentum tensor, to calculate full order corrections to the correlators 
with quantum numbers of nucleon and A-isobar. Expressed in terms of physical nucleon 
and A contributions, these corrections determine the leading low T behavior of the thermal 



mass shifts, which are of order T'* and negative, as in the case of p and oi mesons [ 11 1. The 
mass shifts are rather small numerically, reaching a few MeV at T ~ 100 MeV. This is 
an order of magnitude less than what one would get from the naive scaling, ?n'Ar(A)(^) ~ 
i—{qq)T)^^^, motivated by the results of Ref.|12]. 
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